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To acquire the best path-entangled photon Fock states for robust quantum optical metrology with 
parity detection, we calculate phase information from a lossy interferometer by using twin entangled 
Fock states. We show that (a) when loss is less than 50% twin entangled Fock states with large 
photon number difference give higher visibility while when loss is higher than 50% the ones with less 
photon number difference give higher visibility; (b) twin entangled Fock states with large photon 
number difference give sub-shot-noise limit sensitivity for phase detection in a lossy environment. 
This result provides a reference on what particular path-entangled Fock states are useful for real 
world metrology applications. 
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I. INTRODUCTION 

The application of quantum states of light has long 
been proposed to achieve greater precision, resolution, 
and sensitivity than what is possible classically [TJ |^ . A 
maximally path-entangled state is a superposition of all 
photons in one path with none in the other, and vice 
versa. These states are known as NOON states and were 
introduced to achieve high resolution and high sensitiv- 
ity in metrology and imaging |3]. They are defined 
as |iV :: 0)^.6 = l/%/2(| A^, 0)a,6 + |0, iV)a,fc), where a and 
b indicate the two paths of a two-mode interferometer. 
However, NOON states tend to decohere easily when pho- 
tons are lost from the system. This makes A^OOA^ states 
unusable in real life, where loss is almost always present 
[5H7]. In 2008, Huver et al. proposed a class of gener- 
alized Fock states where decoy photons are introduced 
to the NOON state in both arms of a two-mode inter- 
ferometer [S]. These are called mm' states and they are 

denoted |m :: m')a,b = ^/ V^i\m,m') a,b + \m',m)a,b)- K 
was discovered that mm' states have better metrological 
performance over NOON states in the presence of photon 
loss. In this paper, we locate the best performing m and 
m' under certain fixed loss, where the photon number 
difference (Am = m — m') between the two arms in the 
initial state is fixed. 

In addition to state preparation, achieving super- 
resolution (beating the Rayleigh limit [HI [TO]) and super- 
sensitivity (beating the shot-noise limit [TTj) requires de- 
tection schemes with particular properties. In this paper 
we choose parity detection, which reaches Heisenberg- 
limited sensitivity when combined with lossless NOON 
states [T^HTi] . The parity operator can be expressed as 
H = exp(i7rn) in the number basis or H = exp(i7r(Jo — 
Jz)) in Schwinger notation [T5l - fl7] . which will later be 
discussed in more detail. The parity operator is assigned 
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FIG. 1. A lossy interferometer model which can be used to 
measure the accumulated phase on path b. The phase shift 
operator is given hy U — exp(i</)S^6). 



a parity of +1 if the measured number of photon is even 
and a parity of —1 if odd. In this paper we use the parity 
operator transformed through a beam splitter. Readers 
who are interested in more details about the parity op- 
erator and its application in quantum optical metrology 
may refer to Ref. [H] . 

In this paper we calculate both the visibility and sen- 
sitivity of the phase signal from the interferometer. On 
one hand, the signal from parity detection can be nega- 
tive therefore the ordinary definition of visibility is not 
applicable. We define a "relative visibility" in Section V 
to solve this problem. On the other hand, sensitivity is 
defined by the linear error propagation method as [18j 



AH 



\d{n)/d<i>\' 



(1) 



and AH = y (H^) — (H)^. A promising goal of this work 
is to provide a strategy for choosing the path-entangled 
number state that optimizes either visibility or sensitivity 
for a given loss. 
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II. DENSITY MATRIX 

We start with the classical Mach-Zehnder interferome- 
ter as shown in Fig. [1} where the source and the detector 
are represented by their respective boxes. Similar to the 
approach in Ref. [TH [2^ , the loss in the interferome- 
ter is modeled by adding fictitious beam splitters. Notice 
that it will not change the density matrix if the beam 
splitter is placed before the phase shifter. This can be 
easily demonstrated by removing phase dependence from 
Eq. ([5]) (see below) and applying (j = exp(i(/)6^6) on it 
afterwards. 

The wave function for the mm' input state at stage I 



\m ■■■■ m')a',b' = --j^{\m,m')a',b' + \m' , m) a- .b') ■ (2) 

Without loss of generality, we assume Am = m — m' 
is positive (a mm' state reduces to a NOON state when 
m = N and m' = 0). 

Then the phase shifter introduces a phase shift cj) on 
arm b so that the state at stage II becomes (see Fig. [T]) 



= a\m,m')a',b' + PW ,m) a' ,b' , 



= ^ (e™'^|m, m')a'^b' + e^^m', m),,,^,) (3) 



where a = e™''^/V2 and /3 = e""'^/%/2. We can see that 
because of the different number of photons being phase- 
shifted on arm &, the two paths accumulated different 
phase shifts and thus provide the possibility of interfer- 
ence upon detection. 

The mode transformation by the beam splitter is given 
by Ref. [21] 



b = tlb' + r*,v, 



(4) 



where U = ^/Ji exp{iipi) and — \/Ri exp(j-!/),i) (i — a, b) 
are the complex transmission and reflectance coefficients 
for modes a and 6, where Ti + Ri = 1. 

By tracing out the vacuum modes on both paths, we 
have a density matrix pab that corresponds to the output 
field as 

Pab = Tr„„,„J|*)iii iii(*|] 

m m' 

= ;^^(|apdi|A:,fc')a,ft a,b{k,k'\ 

k=0 k'=0 

+ \P\^d2\k',k)a^b a,b{k\k\ 
m m' 

+ ^ ^ (a/3*d3|Am -I- fc, k')a,b a,b{k, Am -t- k'\ 

k=0 fc'=0 

+ a*l3d4\k', Am + k)a,b a,b{^rn + k', k\), (5) 



where coefficients di{i = 1, 2, 3, 4) are defined as 



dsiKk') = 



If m 



2\Am + fc/ VAm + fc 



^ ^ m'\ ^ ( m'^ ^ 



k' 



dilk, k!) 



If m 



2VAm + fc/ VAm + fc 



i i ^ 

^ m'\ ^ ( m'^ ^ 



fc' 



i(Am+2/c')^,n'-fe' j^5(Am+2fc) j->,„'_fc 



(6) 



An equivalent way to describe the loss process is by using 
the Kraus operators and one may refer to Ref. [22l - f24] and 
references therein. 



III. PARITY OPERATOR 

The calculation of this section is done in Schwinger no- 
tation, so we wish to discuss this representation briefiy. 
Typical four-port two-mode interferometers can be de- 
scribed using Schwinger notation isomorphic to angu- 
lar momentum 15j. The operators are: J^^ — (d^b + 
bU)/2jy = {d^b-b^a)/{2i),J, = {a^ a - b^b)/2, Jq = 
{a'^d + b'^b)/2 and P = Jl + + Jl, which obey 
[Ji.Jj] = iSijkJk and [Jo,Ji] = (i = x,y,z). There- 
fore the common eigenstate of Jq and Jz is the two-mode 
Fock state 



b»2 = \.j + P,j - fJ-)a,b 



(7) 



with eigenvalues j = (Jo) and /i = (Jz)- One uses 
in quantum mechanical angular momentum treatment 
where Jq is not well defined [3S]. However, Jq is use- 
ful in quantum optics because it is directly related to the 
total number of photons in the system. 

In this paper we start with II = (^—lyo-J^ at stage IV 
and transform it back to stage III as Q. The generator 
for the beam splitter transformation is J^, and we have 



7r - - 7r - 

exp(-i- J:r)nexp(i- Ja;) 

exp(i7rJo) exp(i7rjj,). 



(8) 



Following Ref. |17| . the parity operator inside the in- 
terferometer in number basis becomes 



N 



(9) 



fc=0 



where the first summation is over all possible photon loss, 
and N is the total number of photons without loss. It is 
easy to check that = 1. 
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IV. APPLICATION OF PARITY DETECTION 
WITH LOSS AND mm' STATES 

With both the density matrix and the parity operator 
obtained at stage III, it is straightforward to calculate 
the expectation value of the parity operator for an mm' 
state as 



(Q) = Tr(Qp) 

= Ki + if 2 cos Amcj), 

where Ki and K2 are defined as 

m' 

Ki =J2idi{k,k)+d2{k,k)) 

k=0 

(Tim' jym I jym r/m' \ 

7-! / / -1 ^a-^b \ 

X 2Fi[-m,-m ; 1; ^-^j, 

UaKb 

m' 

K2=Y,iMk,k)+d4k,k)) 



(10) 



(11) 



fc=0 



T-}7n Dm rji „ rri - 

^ J^b 



m 
Am 



TaTb 



X2F,{-m',~m';l + Am;j^). (12) 



Here 2-^1(0, b; c; z) 



00 



ia)n{b)n Z" . 



is the ordinary hy- 



•^-^ (c)„ n\ 

n=Q ^ 

pergeometric function [26 . The Pochhammer symbol 
within are defined to be 



{x)r 



1, ifn = 0; 

x{x + 1) . . . (x + n — 1), if n > 0, 



which truncates the infinite summation in the hyperge- 
ometric function at n = m' . Note Eq. (10 1 and (12) 



reduce to the NOON state result if m = and 771' = 0. 

For later calculations and analysis in this paper we use 
loss rate Li = l — T^ {i — a, b) instead of the transmission 
rate Tj following traditional notation in metrology. 



V. VISIBILITY 

We use the parity operator for detection and its expec- 
tation value can be negative in certain regions of param- 
eter space. To quantify the degree of measured phase in- 
formation we need a proper definition of visibility. From 
Eq. ( [T2| ) we can see that Ki decreases and K2 increases 
as the loss rate decreases. Hence, Ki and K2 have a 
range from to 1 and so it is reasonable to define the 
measured signal as 



S = 



K2 
K1+K2 



(13) 
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FIG. 2. (Color online) From left to right, top to bottom. 
Visibility for mm' states with Am — 1 (|1 :: 0) and |3 :: 2)) 
and mm' states with Am = 4 (|4 :: 0) and |6 :: 2)), as a 
function of loss in both arms of two-mode interferometer. La 
and Lb, respectively. Contour lines represent the value of the 
visibility. 



which is always positive. We can then define a visibil- 
ity related to the highest phase information degree (i.e. 
strongest signal) as 



ViLa,Lb 



S{La,Lb 

3(0,0) 



(14) 



where S{0, 0) represents the signal without loss. It is easy 
to see this relative visibility has a value ranging from 
to 1. 



A. Visibility for general cases 

In Fig. [2j we plot visibility as a function of loss rate. 
We see that, for mm' states with a large total number of 
photons, the visibility changes rapidly at high or low loss 
but slowly at mild loss. 

To clearly see the effect of photon number on visibil- 
ity, we assume La ^ Li, = L and plot visibility as a 
function of L for different states in Fig. [3j We observe 
mm' states exhibit a lower visibility than corresponding 
iVOOA^ states for loss rates lower than 50%, and exhibit 
higher visibility for loss greater than 50%. Each row has 
a fixed photon number difference and the total number 
of photons increases from left to right. We can see that 
with increasing total number of photons, the distance 
between the mm' state and the NOON state curves in- 
creases. Each column has fixed to' and the photon num- 
ber difference increases from top to bottom. We can see 
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FIG. 3. (Color online) Visibility for mm' states as a function 
of loss L in both arms of the two-mode interferometer. Each 
row has a fixed photon number difference and the total num- 
ber of photons increases from left to right; each column has 
a fixed m' and photon number difference increases from top 
to bottom. The red dashed lines represent the corresponding 
iVOOiV states with iV = Am. 



that as the photon number difference increases, the dis- 
tance between the mm' state and the A^OOA^ state curves 
decreases. 

Therefore, to obtain the best visibility, under photon 
loss less than 50%, A^OOiV states should be used with as 
many photons as possible, i.e. the bottom right corner of 
the figure; for loss greater than 50%, mm' states should 
be used with as many photons as possible while keeping 
photon number difference to a minimum, i.e. the upper 
right corner of the figure. 

Mathematically, the above results can be explained by 
expanding the visibility of any mm' state around L ^ ^ 
as 



B. Visibility for extreme cases 

For situations with almost no loss, i.e. the loss rate 
L — > 0, the visibility function Eq.(14) can be expanded 
as 



V\ 



L«0 



Am 



OIL 



Am+ll 



(16) 



which explains the behaviors of visibility curves around 
L = for different Am in Fig. [3j Similarly, visibility for 
very lossy situations can be easily expanded as 



V\l^i = 1- 



m 
Am 



(1-L)^™ + 0[(1-L)^ 



(17) 



because of the symmetry of the system. Another exam- 
ple for symmetry is that for 50% loss the visibilities are 
calculated to be exactly one half for all m and m' value, 
which is the consequence of Eq.(15). 



VI. SENSITIVITY 

Another important quantity in quantum optical 
metrology is the precision, or sensitivity, of the phase 
measurement. The Heisenberg limit for any mm' state 
under loss rate La and Lf, should be 1/N while the 

corresponding shot-noise limit is 1/\/n where N — 
[m -t- m'){\ — La/I — L}j/2) is the effective number of 
transmitted photons. Therefore one usually compares 
the performance of different states with the same total 
number of photons. However, in order to keep the same 
resolving power we fix the photon number difference Am 
between two arms of the two-mode interferometer in this 
section. 



A. Sensitivity for general loss 



Sensitivity calculated from Eq. ([T]), ([To]) and (12 1 can 
be expand as 



Ar 



-(L--) + 0[iL--n (15) 



Note that any A^OOA^ state is just the corresponding mm' 
state with m' — 0, therefore it gives the steepest slope 
around L = 1/2 in every case. Physically speaking, this 
result is different from Ref. |8], where NOON states al- 
ways have lower visibility than mm' states under any loss. 
The reason for this discrepancy is that all off-diagonal 
terms of the density matrix are included in Ref. [8| while 
here the parity operator collects only part of the off- 
diagonal terms, making the amplitude of the resultant 
signal smaller. 



_L_ + i!r^csc(Am0)L + O[L^], 
if Am is even; 

^+^-^sec{Am^)L + 0[L% 
if Am is odd. 



(18) 



It is then trivial to show that at the limit of L 0, 
mm' states and corresponding NOON states approach 
minimal phase sensitivity 8(j)rnin = 1/Am at optimal 
phase shifts 



(2rt-l) 7r 
2 A?i 

Am ' 



if Am is even; 

if Am is odd, n — 1,2, . 



(19) 
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FIG. 4. (Color online) Sensitivity for mm! states as a function 
of loss rate L in both arms of a two-mode interferometer. 
Solid blue curve corresponds to |6 :: 0), dashed blue curve to 
l/-y/6{l — L), solid red line to j8 :; 2) and dashed red curve 
to 1/ -^10(1 — L). Parity detection using mm' states is worse 
than shot-noise under around 10% loss. 



FIG. 5. (Color online) Sensitivity of mm' states as a function 
of phase shift (j> from a two-mode interferometer under 35% 
loss. Solid blue curve corresponds to |6 :: 0), dashed blue 
curve to x 0.65, solid red line to |8 :: 2) and dashed 

red curve to 1/ -^10 x 0.65. While the sensitivity of mm' and 
NOON states are worse than respective shot-noise limit, mm' 
is more robust as expected. 



For a mm' state or iVOOA^ state to be able to beat 
the shot-noise limit under parity detection, we therefore 
should have 



Am > Vti 



(20) 



To meet the above criteria, here we choose Am = 6 and 
(j) = , assuming La — Li, for practical purposes. In 
Fig. |4pt can be seen that |6 :: 0) and |8 :: 2) give sub- 
shot-noise performances for loss less than about 10%, and 
|6 :: 0) gives higher sensitivity than |8 :: 2) up to 25% 
loss; for loss greater than 25%, |6 :: 0) is outperformed 
by |8 :: 2) but both are worse than shot- noise. 
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B. Sensitivity for smaller loss 

Often we are more interested in low-loss regions where 
the sensitivity of mm' states and iVOOiV states are com- 
parable to the shot-noise limit. Fig. [5] shows that the 
sensitivity of |6 :: 0) and |8 :: 2) are noticeably worse 
than the respective shot-noise limit under moderate loss 
(35% in this case). Here |8 :: 2) turns out to be more 
robust than |6 :: 0) as predicted in Ref. |S]. 

This robustness, however, does not apply to situations 
where the loss is even smaller. In Fig. |6] we show the 
sensitivity of |6 :: 0) and |8 :: 2) under 5% loss. Here both 
states give higher sensitivity than the shot-noise limit and 
A^OOA^ is the best of all. In contrast, the result in Ref. [8 
shows that, with a certain detection operator, mm' states 
always give better sensitivity than using iVOOA^ states 
and the shot-noise limit no matter how high the loss. 
This discrepancy indicates that the parity operator is not 
the optimal detector that always favors mm' states (see 
conclusion for more discussion.) In addition, more states 
of the form |m :: m — 6) with their highest sensitivity and 
corresponding shot-noise limit are shown in Table [ij 



FIG. 6. (Color online) Sensitivity of mm' states as a function 
of phase shift from a two-mode interferometer under 5% 
loss. Solid blue curve corresponds to |6 :: 0), dashed blue 
curve to l/y/G X 0.95, solid red line to |8 :: 2) and dashed 
red curve to 1/ y/10 x 0.95. While the sensitivity of mm' and 
NOON states are better than their respective shot-noise limits, 
NOON is more robust, which is unexpected. 



TABLE I. Optimal sensitivity of different mm' states 
with Am = 6 and corresponding shot-noise limit 
l/^(m + m/)(l - L) under 5% loss. 



m m/ S(P SNL 

6 0.227 0.419 

8 2 0.266 0.324 

10 4 0.307 0.274 

12 6 0.348 0.242 

14 8 0.387 0.219 
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To conclude, mm' states with large m and small m' 
or NOON states under small loss perform with sub-shot- 
noise limit sensitivity. 



VII. CONCLUSION 

We calculate the visibility and sensitivity of a phase- 
carrying signal by using mm' states with parity detection 
in a lossy environment. In our calculation, we take the 
photon number difference between the two arms of the 
two-mode interferometer Am to be fixed to maintain the 
desired resolving power. Since visibility is not well de- 
fined for signal using parity detection, we use a visibility 
which is measured against signal without loss. To have 
high visibility, one should use NOON states with large N 
(when loss is low) and mm' states with large m + m' but 
small Am (for high loss). Considering only sensitivity, 
our calculation shows that NOON states with large N or 
mm' states with large 5m under low loss are capable of 
performing with sub-shot-noise limit precision. 

It is worthwhile to mention two points. First, modeling 
loss with a single fictitious beam splitter is sufficient for 
practical purposes. If one reverses the order of the phase 



shifter and the beam splitter in Fig. [Tj the same density 
matrix will be obtained as proved in Ref . [23l [24] . This 
also means our model is equivalent to a continuous loss 
model. Second, Ref. [5] uses a detection operator that is 
carefully chosen so that it sums up all off-diagonal terms 
of the density matrix and provides sub-shot-noise sensi- 
tivity. Meanwhile the parity operator collects some of the 
diagonal and off-diagonal terms. The inclusion of diago- 
nal terms may reduce the signal size and therefore visi- 
bility or sensitivity of phase information. However, such 
an operator in Ref. [5] is yet to be produced in a lab set- 
ting. On contrary, a lot of effort has been made to realize 
parity measurements. A straightforward parity measure- 
ment relies on high precision photon number-resolving 
detection at single-photon level, which has been demon- 
strated experimentally in near-infrared region ^27j . Alter- 
native parity detection setups without number-resolving 
detectors have been proposed as well [28l ES] . 
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